Bethe-Salpeter equations for the collective modes of the t — U — V — J model with 

d-wave pairing 
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The Bethe-Salpeter equations for the collective modes of a t-U-V-J model are used to analyze 
the resonance peak observed at Q = (tt, 7r) in neutron scattering experiments on the cuprates. 
We assume that the resonance emerges due to the mixing between the spin channel and 19 other 
channels. We have calculated the energy of the lowest mode of the extended Hubbard model ( J = 0) 
vs the on-site repulsive interaction U, as well as the UJ lines in the interaction parameter space 
which are consistent with the ARPES data and reproduces the resonance peak at 40 meV in Bi2212 
compound. We find that the resonance is predominantly a spin exciton. 

PACS numbers: 71.10.Ca, 74.20. Fg, 74.25.Ha 
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Introduction. It is widely accepted that: (i) the 
angle-resolved photoemission spectroscopy (ARPES) 
data produce evidences for the opening of a d-wave pair- 
ing gap in cuprates compounds described at low energies 
and temperatures by a BCS theory, and (ii) the basic 
pairing mechanism arises from the antiferromagnetic 
exchange correlations, but the charge fluctuations associ- 
ated with double occupancy of a site also play an essential 
role in doped systems. The simplest model that is con- 
sistent with the last statements is the t-U-V-J model. 
In the case of d-pairing the gap function is Ak = Ad-^/2, 
where A is the maximum value of the energy gap and 
= (cosk x — cos ky) (lattice constant a = 1). The 

BCS gap equation is 1 



2 J — 7T J — 7T 



where = 2V + 3J/2, E(k) = + A 2 . The 
mean-field electron energy £k has a tight-binding form 
£4c = t\ (cos k x + cos k y ) /2+t2 cos k x cos k y -\-t^,{cos 2k x + 
cos2fc J/ )/2 + ti (cos 2k x cos k y + cos2fcj,cosfc a; )/2 + 
t$ cos 2k x cos 2k y — /i obtained by fitting the ARPES 
data with a chemical potential fi and hopping amplitudes 
ti for first to fifth nearest neighbors on a square lattice. 
A, tx, t 5 and \i should all be thought of as an effective 
set of parameters, while V$ has to be determined by 
the gap equation. For Bi2212 compound, there are two 
possible sets of parameters with all tight-binding basis 
functions involved (see Table 1 in Ref. [l|). Assuming 



35 meV, we obtain V 



(i) 



115.2 meV with set 1, 
and V^ 2) = 87.9 meV with set 2. Hao and Chubukov 2 
have used another set of parameters (we shall call it 
H&C) for Bi2212 compound with a doping concentration 
x = 0.12: tx = -At,t 2 = 1.2t, t = 0.433 eV, /z = -0.94*, 
A = 35 meV and V$ = 0.6t. The parameters U, V 
and J should be adjusted in such a way that the 
sharp collective mode which appears at wave vector 
Q = (tt, it) in inelastic neutron-scattering resonance 
(INSR) studies- occurs at energy which corresponds 
to the lowest collective mode of the corresponding 
Hamiltonian. In RPA the resonance is determined by 
the pole of the spin correlation function, which in the 
case of J = (the phase diagram at half filling shows an 



"island" in U-V space where d-wave pairing exists^) is: 

XsM = XoolQo^)/! 1 + UXooiQoi")}' where the bare 
spin correlation function^ is Xoa = -^77 (^77 ^ s defined 
later in the text). Using the H&C set of parameters^ 
and a resonance energy of 40 meV, we calculate the 
RPA value of U of about 1.16 eV. Sets 1 and 2 provide 
U {1) = 0.533 eV and U {2) = 0.418 eV, respectively. The 
coupling of the spin channel with other channels should 
change the RPA results for U. For example, we have two 
7T channels^ with bare it susceptibilities Xxi = ^ii an< ^ 
X22 = ^77; respectively. The susceptibilities I^~,Jj~,Jff 
represent the mixing of the spin and two it channels. 
Thus, the coupling of the spin and two tt channels (a 
three-channel response-function theory) leads in the 
generalized random phase approximation (GRPA) to a 
set of three coupled equations and the value of U is 
reduced from 1.16 eV to 0.974 eV. When the extended 
spin channel is added to the previous three channels, 
we have a set of four coupled equations (a four-channel 
theory), and according to Ref. 7] U ~ 300 meV is 
required in the case when V^,=0.260 eV and J = 0. 

In what follows, the energy of the resonance is ob- 
tained from the solution of 20 coupled Bethe-Salpeter 
(BS) equations for the collective modes in GRPA, i.e. 
the resonance emerges due to the mixing between the 
spin channel and other 19 channels. In our approach 
the INSR energy solves detlx^ 1 — V\ — 0, where the 
mean-field response function x an d the interaction V are 
20 x 20 matrices. The secular determinant can be rewrit- 

A B 



ten as det\x 1 — V\ — det 



B C 



= det\C\det\A - 



BC~ X B T \. In the case of the four-channel response- 
function theory^' 8 A is a 4 x 4 matrix while the mix- 
ing with the other 16 channels is represented by a 4 x 4 
matrix BC~ 1 B T . We emphasize that none of the previ- 
ous theoretical interpretations of the INSR feature at Q 
have accounted properly for the mixing term BC~ 1 B T . 

t-U-V-J model. The Hamiltonian of the t-U-V-J 
model consists of t and U terms representing the hopping 
of electrons between sites of the lattice and their on-site 
repulsive interaction, as well as the spin-independent at- 
tractive interaction V and the spin- dependent antiferro- 



2 




U [eV] 



FIG. 1: The energy of the resonance obtained from the BS 
equations when J = 0. The curves are plotted using param- 
eters given in Table 1 in Ref. [U] (set 1 and set 2), and the 
Hao and Chubukov parameters (curves H&C). The puncture 
curve represents the three-channel energy (Fig. 4 in Ref. 



magnetic interaction J: 



— V ni t<T rij t(T > + J ^ S».Sj. 



(1) 



<i,j>crcr' 



<hJ> 



Here, the Fermi operator ^>j CT (V'i.er) creates (destroys) 
a fermion on the lattice site i with spin projection 
(T =t,i along a specified direction and n%, a = ip\ a ipi t<T 
is the density operator on site i with a position vec- 
tor Yi. The symbol X)<ij> means sum over nearest- 
neighbor sites, tij is the single electron hopping in- 
tegral. The antiferromagnetic spin-dependent interac- 

'<i,j> 



tion JJ2<i 7> 

= Ji + J2 consists of two terms: 

Ji 

and J2 = J 2 



lE<y>Mj,t 



2 <0<i,j> 

It is useful to 



introduce four-component Nambu 
fermion fields tp(y) = ^(^{(^T^V^y)) and 

4>{x) = (^p^(x)ip^(x)')p^(x)^pi(x) S j , where x and y are 

composite variables and the field operators obey anticom- 
mutation relations. The "hat" symbol over any quantity 
O means that this quantity is a matrix. 

The interaction part of the extended Hubbard Hamil- 
tonian is quartic in the Grassmann fermion fields 
so the functional integrals cannot be evaluated ex- 
actly. However, we can transform the quartic 
terms to a quadratic form by applying the Hubbard- 
Stratonovich transformation for the electron operators: 9 

J DAe ^A a (z)D^-\zy)A fi (z)+^(y)fi°Hy-,x\z)ip(x)A a (z)] = 



last equation is used to define the 4x4 matrices D^l and where G(k; w) is the BCS Green's functionJ^iii The di- 



The 



FIG. 2: Line in U, J parameter space which reproduce the 
INSR energy of 0.04 eV. Note that V = - 3J/4 where 

Vip = 0.6i = 259.8 meV is calculated from the gap equation 
by using the set of parameters given in Ref. 0. 



Ta* 1 (cx,/3 — 1,2,3,4). Their Fourier transforms, written 
in terms of the Pauli <7j, Dirac 7 and alph a 10 i 1:L matrices, 



are as follows: = 



D l 
D 2 



f<°> = ( 7 °±a,)/2 



and = (a x ± ia y )/2, where a. 



<J l 
<7,.<j,-<7„ 



Bi = ( J(k) - V(k)) a-o + (U - J(k) - V(k)) a x and 
D-2 = 2 J(k)<j x . For a square lattice and nearest-neighbor 
interactions V(k) = 4V(cos(k x ) + cos(k y )) and J(k) = 
J(cos(fc x ) +cos(fc y )). Now, we can establish a one-to-one 
correspondence between the system under consideration 
and a system which consists of a four-component boson 

field A a (z) interacting with fermion fields ip(y) and ?p(x). 

(e) (A) 

The action of the model system is S — S + S Q + 

%)GM-i(y;x)f(x), S { 
A) _ 7T,f».\r(°) 



S^-^ where 
1 a /^nl )- 1 



(A) 



2 A a {z)D^>-\z,z')Ap{z') and S^~ A ) = ^ (y, x 

z)ip{x)A a {z). Here, we have used composite variables 
x,y,z — {ri,u}, where is a lattice site vector, and 
variable u range from 0to/3 = l/fcsT'(T and fee are the 
temperature and the Boltzmann constant). We set h = 1 
and we use the summation-integration convention: that 
repeated variables are summed up or integrated over. 

Following the same steps as in Refs. [H 13 1 , we can de- 
rive a set of sixteen BS equations for the collective mode 
w(Q) and BS amplitudes *£ in2 (k) (m,n 2 = 1,2,3,4). 
Their matrix representation at zero temperature is : 



rW*Q(q)F^G(k; fi) + D { a °>(Q)G(k + Q; + u ) 
fi )G(k;^)Tr[f( 0) $Q(q)]}, 



(2) 
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FIG. 3: Line in U, J parameter space which reproduce the 
INSR energy of 0.04 eV. The V value is V = - 3J/4 

where = 115.2 meV and V^ 2) = 87.9 meV are calculated 
by using two sets of parameters given in Ref. 0]. 



rect D^l (k-q)rL 0) f f and exchange 5gj (Q)f f in- 
teractions mix all sixteen BS amplitudes. We can greatly 
simplify Eqs. @ using the fact that in the RPA the sus- 
ceptibilities at Q are convolutions of two single-particle 
Green's functions G, and the equation for the collective 
mode in the RPA is: V)X2° )_ V) - C7ia(w) = 0, 
where the susceptibilities Xi an d X2 originate from 
(U,Ji) and J2 interactions, respectively. The term G12 
mixes the J\ and J 2 interactions, but it is proportional to 
convolutions which involve the anomalous Green's func- 
tions G13 and G24. The two Green's functions appears 
in the case of spin triplet pairing states where the order 
parameter A a!| g(k) is a 2 x 2 matrix. For a singlet super- 
conductivity and d-wave pairing A Q/ 3(k) = i(a l/ ) Q( gA(k), 
Ci2(w) = 0, and the equation for collective modes be- 
comes [1 + (U + 4J).%] [1 + 4J.%] = 0, i.e. Ji and 
J 2 terms contribute separately to the collective modes. 
Thus, we shall neglect all contributions due to the J 2 
term in Eqs. ([2|). In this approximation we have a set of 
four equations, which can be further simplified to a set 
of two equations in the same manner as in Refs. [l2lfl~3| : 



w(Q) - e(k, Q)]G+(k, Q) - ^ £ [ 7 k, Q 7 q ,Q + /k, Q / q ,q] G+(q, Q) - ^ ]T [7k,q7 q ,q - ik.qiq.q] G~(q, Q) 

q q 

- E t y ( k + JO* - q)] [7k, Q 7q,q + Kq^,q] G + (q, Q) 

q 

~ 2N E ^ k ~ q ^ ~ J( ' k ~ q ^ ^ kQ ^ qQ + m k,Q"i q:Q ] G+(q, Q) 
q 

+ ^ E t y ( k - 1) + J ( k - <J)] [7k,Q7q,Q - «k,Q« q , Q ] G-(q, Q) 
q 

+ 2A^ E ~ q ) ~ J ( k ~ ^k,Q7q,Q - TO k,q m q,q] G"(q, Q) 

q 

- U ~2N Q) E ^.Q7q,Q (G + (q, Q) - G- (q, Q)) + - ~ 2 ^ (Q) ]T m k , Q m q , Q [G+ (q, Q) + G~ (q, Q)] , 



U 



(3) 



[w(Q) + e(k, Q)]G (k, Q) = ^ [7k, Q 7 q ,q + k.q/q.q] G (q, Q) + ^ £ [^,Q7q,Q - ^,q^,q] G+(q, Q) 

q q 

+ ^ E t y ( k - q ) + J ( k - q )] [7k,Q7q,Q + Wq,Q] G~(q, Q) 
q 

+ 2A? E ( k ~ q ) ~ J ( k ~ q )l Rk,Q7q,Q + ^k.QTOq.q] G~ (q, Q) 
q 

- 2^ E t y ( k - q ) + J ( k - q )l [7k,Q7q,Q - *k, Q Zq,Q] G+(q, Q) 

q 

- J^Yl [ y ( k_ q ) _ J ( k ~ q)l [7k,Q7q,Q - ^k,QTO q ,q] G+(q,Q) 

q 

- U ~2N Q) E 7^Q7q,Q (G + (q, Q) - G~ (q, Q)) - U "^ (Q) £ m k , Q m q , Q [G+ (q, Q) + G~~ (q, Q)] . 



(4) 



4 



Here e(k, Q) = P(k-I-Q) +P(k), and we use the same form factors as in Ref . [12| : 7k, q = Uk"k+Q + Vk^k+Q, ^k.Q = 
Wk"k+Q-w k Wk+Q, 7k,Q = "k^k+Q-Wk+Q-Uk, and m k , Q = u k w k+ Q+u k+Q i; k where u 2 . = l-v£ = [1 + e(k)/£'(k)] /2. 

It is worth mentioning that in the case of an extended Hubbard model (J = 0), Eqs. (j3|) and ((4|) are the exact BS 
equations in the GRPA. They are in accordance with the Goldstone theorem which says that the gauge invariance 
is restored by the existence of the Goldstone mode whose energy approaches zero at Q = 0. The last statement 
corresponds to the so-called trivial solution of the BS equations: G + (k, Q = 0) = — G~(k, Q = 0) = A k /2P(k), and 
the gap equation^ A k = J2ql~U + ^C 4 — A q /2P(q) is recovered from our BS equations. 

The Fourier transforms of V and J interactions are separable, i.e. V(k — q) = 2V^A k Aq and J(k — q) = JA k Aq/2, 

and therefore, Eqs. ([3]) and Q can be solved analytically. Here Ak = (s k , <i k , ss k , sd^) is an 1 x 4 matrix, and 
we have used the following notations: s k = cos(fca;) + cos(fcj,), c? k = cos(k x ) — cos(k y ), ss k = sin(A; x ) + sin(k y ) and 
cdk = sin(fc :c ) — sm(ky). Thus, we obtain a set of 20 coupled linear homogeneous equations for the dispersion of the 
collective excitations. The existence of a non-trivial solution requires that the secular determinant c?ei||x _1 — V\\ 

_ ( P Q 

is equal to zero, where the bare mcan-ficld-quasiparticle response function v = „ T „ 

y V R 

V = diag{U, U, - (U + 4 J), U+ 16V, -(2V+J/2), -(2V + J/2), -(2V- J/2), -(2F- J/2)) are 20 x 20 matrices. 
P and Q are 4x4 and 4 x 16 blocks, respectively, while R is 16 x 16 block (in what follows i, j = 1, 2, 3, 4): 



and the interaction 
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The quantities i" 0i (, = P a .b(e(k, Q)) and J a .b 
&(k, Q) = ?k, q, m k ,Q, 7k,Q or 7u,q): 



and the 4 x 4 matrices FJ b = F^ b (e(k,Q)) and J^ ( 



F a<b (uj), the 1 x 4 matrices P^ = i* 6 (e(k,Q)) and j; 6 = j* 6 (w), 
F % J b {uj) are defined as follows (the quantities a(k, Q) and 



P a ,b(z) 



1 



\ " 



k 



a;a( k, Q)6(k, Q) 
-e 2 (k,Q) 



,FUx) 



N ^ 

k 



zo(k, Q)6(k, Q)A k 
2 -e 2 (k,Q) 

I 



TV 4-^ 



za(k,Q)6(k,Q) 

CJ 2 -£ 2 (k,Q) 



(AjA k ) i 



The elements of P, Q and P blocks are convolutions of 
conventional two normal GG, two anomalous FF Green's 
functions or FG terms. At the high-symmetry wave vec- 
tor Q , P a b and P a b with i = 3, 4 involve sine functions, 
and therefore, all vanish. I 2 b and J 2 6 also vanish because 
e(k, Qo) is symmetric with respect to exchange k x O fey. 
Similarly, the non-diagonal elements of I^ b and J^ b with 
i 7^ j all vanish. Thus, blocks P and Q, each has 10 differ- 
ent non-zero elements, while R has 40 non-zero elements. 
In other words, the ui dependence of % (or x _1 ) comes 
from these 60 non-zero elements. It is worth mentioning 
that within the four-channel theory^ the collective mode 
energy has been calculated by using a 4 x 4 symmet- 
ric matrix \ which has only 6 non-zero elements at Qo: 

Xll = Pr?>X22 = ^ OTm ,X33 = Pj ! 7 \X44 = I?? -,Xl2 = J TO7 

and £34 = Jf^ (the other 4 elements £13 = ^ 77 ,Xi4 = 
J/ 2 75 X23 = Jm~ ( and £24 = I™i vanish). 

In Fig. 1 we present the results of our calculations 
of the lowest collective mode of the extended Hubbard 
model (J = 0) using 49 x 49 k points in the Brillouin 
zone and three possible sets of parameters: sets 1 and 2 
include all tight-binding basis functions (see Table 1 in 



Ref. [![), while the third set (H&C) is used by Hao and 
Chubukov* 2 - As can be seen in Fig. 1, BS equations pro- 
vide energies which are significantly different from those 
obtained according to the three-channel theory (see Fig. 
4 in Ref. 0). In Fig. 2 and Fig. 3 we present the results 
of our calculations of the lines in U, J parameter space 
which reproduce the INSR energy of 40 meV using all 
twenty channels. We see that the RPA spin correlation 
function and the BS equations in GRPA, both provide 
very similar results for U at point J = 0. This indicates 
that the resonance remains predominantly a spin exciton. 

In summary, we have derived a set of four coupled BS 
equations for the collective modes of the t — U — V — J 
model including the J\ part of the antiferromagnetic in- 
teraction. These equations have been used to analyze 
the resonance peak in Bi2212. It is interesting to note 
that the trivial solution of the BS equations © and (|3] 
leads to an equation similar to the gap equation but with 
= 2V+J/2 instead of = 2V+3J/2. The Goldstone 
mode, which is expected on physical grounds as the sym- 
metry is spontaneously broken by the condensate, does 
exist as a trivial solution of the sixteen BS equations. 
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